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Abstract 

We carry out simulations of the anisotropic uniformly frustrated 3D XY 
model, as a model for vortex line fluctuations in high T c superconductors. 
We compute the phase diagram as a function of temperature and anisotropy, 
for a fixed applied magnetic field B. We find that superconducting coherence 
parallel to B persists into the vortex line liquid state, vanishing at a T cz above 
the melting T m . Both T cz and T m are found in general to lie well below the 
cross-over T C 2 from the vortex line liquid to the normal state. 

I. INTRODUCTION 

From a phenomenological point of view, "high T c " superconductors are believed to differ 
from conventional type II superconductors primarily because of the dramatically enhanced 
importance of thermal fluctuations]]]! In an applied magnetic field H, such thermal fluctua- 
tions are believed to melt the Abrikosov vortex line lattice at a temperature T m well below 
the mean field T C 2 which marks the onset of strong diamagnetism.H@ In between T m and 
T c2 is a new vortex line liquid state. Experimentally, this picture has been supported by 
the observation that, in high T c materials, the onset of reversible diamagnetism occurs at a 
temperature well above that where resistance vanishes;! the separation between these tem- 
peratures increases with increasing H. According to this picture, the onset of diamagnetism 
at T c2 is associated with a growth in local superconducting correlations, giving rise on short 
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length scales to a finite superconducting wavefunction ip(v) in terms of which vortex lines can 
be defined. This T c2 marks a strong cross-over region, rather than a sharp thermodynamic 
transition. In the resulting vortex line liquid, free diffusion of vortex lines gives rise to "flux 
flow" electrical resistance. The vanishing of resistance only occurs at a lower temperature 
when the line liquid freezes into a lattice or glass. 

To investigate the effect of thermal fluctuations on phase transitions in type II supercon- 
ductors, within a numerical simulation, one of us and Y.-H. Li have previously introducedi'B 
the three dimensional (3D) uniformly frustrated XY model. Simulations of this model in the 
isotropic coupling limit, at low vortex line density, gave the surprising result that supercon- 
ducting coherence parallel to the applied magnetic field appeared to persist above the vortex 
line lattice melting temperature, into the vortex line liquid phase.i The goal of the present 
work is to extend these simulations to a model with uniaxial anisotropic couplings, so as 
to better model the layered structure of the high T c materials. We will consider only the 
case where the applied magnetic field is parallel to the anisotropy axis z. We will map out 
the phase diagram in the anisotropy-temperature plane, looking for the presence of parallel 
coherence in the vortex liquid phase and dimensional cross-over as anisotropy increases. The 
rest of this paper is organized as follows. In Section |TJ we describe our model, its limits of 



validity, and the specific parameters of our simulations. In Section |TTT| we give our numer- 
ical results, mapping out the phase diagram, and characterizing the nature of vortex line 



fluctuations. In Section IV we discuss our results and present our conclusions. 



II. MODEL 



Our model arises from the familiar Ginzburg-Landau (GL) free energy functional with 
the following additional approximations: (i) the continuum is discretized to a periodic grid 
of sites, (ii) the amplitude of the superconducting wavefunction is taken to be constant 
outside of a vortex core (London approximation), and (Hi) the internal magnetic field B is 
taken to be a spatially uniform fixed constant. These yield the Hamiltonian,@ 
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n[9i\ = - Jij cos(^ - Q 3 - Aij) , (1) 

(i3) 

where 9i is the phase angle of the wavefunction at site i of the discrete grid, the sum is 
over nearest neighbor bonds (ij), Aij = (27r/0o) // A • di is the integral of the magnetic 
vector potential across the bond (0o = hc/2e is the flux quantum), and = J z or Jj_ is 
the superconducting coupling in the direction of the bond. We consider here a cubic grid 
of sites. If we identify the discrete spacing along z as the distance d between copper-oxide 
planes, and the discrete spacing in the xy plane as representing the short length cutoff for 
a vortex core, given by the coherence length we have 

x " 167r 3 Ai ' * IQ^Xld ' K > 

where and X z are the magnetic penetration lengths within and normal to the copper-oxide 
planes, respectively. We define an anisotropy parameter 77 as 

J± X z d 

Note that if the coherence length along z is larger than the interplanar spacing, £ 2 > d, 
then one should replace d by £ 2 in Eqs.(Q-0) above. In this case, since the GL free energy 
with anisotropic masses givesQ £ 2 /£j_ = \±/X z , we have r] = 1 and hence isotropic couplings. 
This isotropic model we have investigated previously.0 In the present paper we extend these 
studies to the anisotropic case 77 > 1. 

In the model of Eq.(|l|) spatial variations and fluctuations in the internal magnetic field 
are ignored, with V x A = Bz a uniform constant. This should be valid provided B is so 
large that the magnetic fields associated with each vortex line strongly overlap, i.e. 



y/B « \ ± , (4) 

where a v is the spacing between vortex lines. However B should still be small enough that 
a v ^> £_l (so details of the vortex cores are not important). The anisotropy must also be 
such that 
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d < \ 2 J\ Z , 



or equivalently, 



(5) 



so that the Josephson coupling between the planes dominates over the magnetic coupling.BB 
Formally, our model corresponds to the limit of Aj_/£_i_ — > oo, keeping J± constant. For our 
simulations, we take the A^ as fixed constants, chosen to give a particular fractional density, 



of vortex lines penetrating the xy plane. 

Using the model of Eq.(|l|), which is in terms of the phase angles 9i, we will also study 
vortex line fluctuations. To locate a vortex line, we compute the phase angle difference 
[8i — 6j — Aij] across each bond, restricting this angle to the interval (— tt, 7r]. The circulation 
of these angle differences around any plaquette a must then be 27i(n a — f a ), where f a = 
or / depending on the orientation of the plaquette, and a non-zero integer value of n a 
indicates the presence of a vortex line piercing the plaquette. Computing the vorticity of 
each plaquette in this fashion, we can then trace out the paths of the vortex lines. 

To model a particular material, we would like to map out the phase diagram as a function 
of T and magnetic field B, for a fixed value of anisotropy rj. However, due to commensu- 
rability difficulties between the triangular vortex lattice preferred in a continuum and the 
discrete sites permitted by our numerical grid, different vortex line densities would form 
lattice structures of differing symmetry in the ground state. Since we are computationally 
limited to a fairly coarse grid, this would make direct comparison of systems with differ- 
ent B difficult. We therefore choose to map out the phase diagram as a function of T and 
anisotropy rj, for fixed B. We can see however, using dimensional arguments, that increasing 
77 at fixed B, is similar to increasing B at fixed rj. If we measure any transition temperature 
T c in units of J±, then the dimensionless T c / Jj_ can only depend on the other dimensionless 
parameters of the Hamiltonian (|IJ), the anisotropy rj = X z d/X±C,±, and the vortex line den- 
sity / = B£j_/(fi Q = (t;±/a v ) 2 . Since our London approximation ignores details of the vortex 
cores, if we consider the continuum limit of our model, a v 3> £±, we expect that T c / J± 
should be at most weakly dependentS on the vortex core radius The only combination 
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/ = Be ± /<Po , 



(6) 



of 7] and / that is independent of £j_ is rf f . Thus, the dominant dependence of T c /J± on 77 
and / can only be through some function of rff ~ ifB. 

We can further argue how transition temperatures should depend on the quantity rff . 
In the limit of extreme anisotropy, 77 — > 00, we have completely decoupled planes, and 
the transition temperature should be independent of 77; thus we expect T c ~ J±. In the 
limit of a nearly isotropic system, 77 ~ 1, we expect that T c should be independent of the 
spacing between planes d; thus we expect T c ~ J_i_/r]\/f = (<pl/16Tr 3 X 2 L )(Xj_/X z )((j)o/B) 1 ^ 2 . 
These are in fact the predictions for the melting temperature based on Lindemann criterion 
calculations .B UB 

The cross-over from small to large 77, where the discreteness of the layering along z 
becomes important and one approaches the two dimensional limit, can be estimated by the 
criterion r) 2 T f ~ 1, or, using / = B^j_/4> = (£j_/a v ) 2 , as 

Vcr = a v /£_ L (7) 

Using an effective elastic medium approximation to describe vortex line fluctuations in the 
line lattice, one can show01 that for 77 < r] CI , the dominant wavenumber q z of fluctuations at 
melting satisfies the condition d < ir/q z , and hence the layering of the material is averaged 
over. For 77 > r] CT however, the dominant wavenumber is at d = n/q z , and layering is 
important. Some theoretical models!'01ll have predicted that i] cr (or equivalently B CT = 
4>q\\/ \ 2 z d 2 ) will mark a dramatic change in behavior, reflecting a three dimensional to two 
dimensional cross-over. Looking for any such cross-over behavior at 7] CT will be one of the 
goals of this work. 

Our simulations consist of standard Metropolis Monte Carlo simulations of the Hamilto- 
nian ([[]), using periodic boundary conditions in all directions, on grid sizes L\ x L z . We use 
a magnetic field B which yields a fractional density of vortex lines / = 1/15. The ground 
state vortex lattice, shown in Fig. ^ is a nearly triangular vortex line lattice with sides of 
length -\/l8 x v^18 x y/Tf in units of To map out the 77 — T phase diagram, we have done 
simulations varying T at different fixed values of 77 on lattices of size 15 3 . We have also car- 
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ried out simulations of larger system sizes for the specific cases of rf = 10 < r)% r = 1// = 15, 
and ?7 2 = 50 > 77^.. Our runs are typically 10, 000 sweeps through the grid to equilibrate, 
followed by 128,000 sweeps to compute averages. These simulations are about 9 times longer 
than in our previous work.0 Errors are estimated by a standard data blocking procedure. 

III. NUMERICAL RESULTS 
A. Phase Diagram 

To test for superconducting coherence, we compute the helicity moduli Tj_(T) and T X (T) 
which measure the stiffness with respect to applying a net gradient ("twist") in the phase 
angle of the wavefunction along directions perpendicular and parallel to the applied magnetic 
field! The helicity modulus in direction /j is given by the phase angle correlation 



T M (T) = 




where is the unit vector from site i to j. When T M is positive, the system can carry a 
supercurrent, and so possesses superconducting coherence in direction ft. When T M vanishes, 
superconducting coherence is lost. 

To determine the vortex line lattice melting temperature, we compute the density-density 
correlation function of vortices within the same plane, 

S(k ± ) = ±- £ ^< r ^(n z (v ± ,z)n z (v' ± ,z)) , (9) 

r ±i r X> z 

where n z (r±,z) is the vorticity at site r_i_ in the xy plane at height z (henceforth, we will 
refer to the vortices in the xy planes as the "pancake" vortices). Below melting, we expect 
to see a periodic array of sharp Bragg peaks in the plane. Above melting, we expect to 
see the broad circular rings characteristic of a liquid. 
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We also measure the specific heat per site of the system, C, using the usual energy 
fluctuation formula. A peak in C locates the temperature at which, upon cooling, there is a 
dramatic freezing out of thermal fluctuations and the system looses the bulk of its entropy. 
We will take the location of a high temperature peak (above any phase transitions) in C as 
indicating the cross-over temperature T c2 where the superconducting wavefunction develops 
on small length scales, vortex lines become well defined objects, and one has the onset of 
strong diamagnetism. 

In Fig. H] we show our results for Yj_ and T z for the case rf = 10. We see that Yj_ vanishes 
at a T C _|_ significantly lower than the T cz where T 2 vanishes. We show data for heating and 
cooling, for three different grid sizes, 15 3 , 30 3 , and 15 2 x 120. Comparing heating and cooling, 
we see no appreciable hysteresis for T 2 . Hysteresis in Yj_ appears only for the 30 3 system, 
where we failed to cool back into a lattice. There are no obvious shifts in T cz or T c ± due to 
finite size effects as L± and L z are varied. We determine our estimates T c ±/J± ~ 0.36 and 
T cz /Jj_ — 0.58 by visually extrapolating the curves to zero from the inflection point that 
marks the onset of the high temperature tails. We have found that the size of these tails 
tends to decrease with increasing simulation time, as well as with system size. 

It is important to note that the finite T c j_ in our model is strictly an artifact of the 
discretizing grid, which acts like an effective periodic pinning potential for the vortex lines. 
In a continuum model, one would find Yj_ = at all temperaturesjl! as the vortex line 
lattice is free to slide as a whole, giving "flux flow resistance." A discretizing grid removes 
this translational symmetry, resulting in a commensurately pinned vortex line lattice at low 
temperatures, with > 0. For a high density of vortex lines, it is likely that the vortex 
lattice remains commensurately pinned until it melts. In such a case one expects T m = T c ±. 
However recent simulations,&"0 with a more dilute vortex line density than studied here, 
have claimed evidence for a depinning T c j_ which is lower than T m , with the intermediate 
phase a floating vortex line lattice. It is thus important to determine the melting T m of our 
vortex lattice independently from our measurement of Tj_. 

In Fig. |^ we show intensity plots at various temperatures of S^kjJ in the kj_ plane, for 
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the 30 3 system upon heating. Looking at when the Bragg peaks disappear, we estimate the 
melting temperature to be T/Jj_ ~ 0.43, somewhat higher than T c j_/Jj_ ~ 0.36. To try to 
quantify the location of the melting transition, we now look at the heights of the Bragg peaks 
at the reciprocal lattice vectors. We denote by {Ki} the six, almost equal, smallest non- 
zero reciprocal lattice vectors. Let {K^} be the six vectors obtained by reflecting the {K^} 
through the x axis. Since the vortex line lattice breaks this reflection symmetry of the square 
discretizing grid, we will have S'(Ki) > S^K'J for the lattice phase. However, once the lattice 
has melted, the reflection symmetry of the grid should be restored. We can therefore define 
as an order parameter of the melting transition AS'(Ki) = S'(Ki) — S(K[). Normalizing 
by So = S(K = 0), and averaging over the six {K^}, we plot in Fig.[| AS'(K 1 )/S'o versus 
T, for the three system sizes, 15 3 , 30 3 , and 15 2 x 120. AS'(Ki)/ Sq decreases linearly over a 
large intermediate range of T. From the 30 3 system we estimate T m / J± ~ 0.44. Note that 
there is a greater finite size effect and more hysteresis for AS(Ki)/Sq than there is in T± jZ . 
The estimate for T m tends to decrease as L± increases. Our result T c ± < T m suggests the 
presence of a floating vortex line lattice. However it remains possible that T c ± and T m will 
merge as the system size increases, due either to the finite size dependence observed in T m , 
or to the possibility that T c ± actually lies farther out in the high temperature tail of Yj_ 
than we have estimated. 

Finally, in Fig. [5] we show the specific heat C. The high temperature peak in C at 
T/J± ~ 1.0 we identify with the cross-over T c2 , which is thus seen to lie well above T cJ _, T m , 
and T cz . A suggestion of a smaller peak is seen at the lower temperature T c ^. 

In Fig. ^| we show Tj_ and T z for the case rf = 50, for system sizes 15 3 and 30 3 . Here 
the data has considerably more scatter than in Fig. ^] (in general, we found it increasingly 
difficult to achieve good equilibration as rj increased). Nevertheless, there again appears to 
be two distinct transitions, with T c ±/J± ~ 0.19 < T cz /J± ~ 0.24. Intensity plots of the 
structure function S'(k^) are shown in Fig.|7|, and the peak height differences AS'(K 1 )/S'o in 
Fig.^|. These suggest a melting T m / Jj_ ~ 0.21. In Fig. |S] we show the specific heat C. The 
high temperature peak at T/J± ~ 1.0 is again associated with T c2 . However, comparing 
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with Fig. [|, there is now a more clearly defined smaller peak at T c ±. 

Carrying out simulations at other values of 77 on a 15 3 grid, we show in Fig. ITU] the resulting 



phase diagram in the rj — T plane. The T cz line denotes the loss of phase coherence parallel 
to the applied magnetic field, as measured by the vanishing of T z . The T m line denotes the 
melting of the vortex line lattice, as measured by the vanishing of AS'(K 1 )/S'o. The T c j_ line 
denotes the depinning of the vortex line lattice from the discretizing grid, as measured by 
the vanishing of Tj_. We see that T m coincides with T c ± in the 77 ~ 1 and 77 ^> r? cr limits, but 
is somewhat greater than T c ± in the vicinity of the cross-over anisotropy r\ cr = a v /^± = \flb. 
Between T m and T cz we have a vortex line liquid which retains superconducting coherence 
in the direction parallel to the applied magnetic field. The dashed line locates the high 
temperature peak of the specific heat, and marks the cross-over from the the vortex line 
liquid to the normal metal. The dotted lines labeled £ c = n will be explained at the end of 
the following section. 

Thus, for T c2 < T we have the resistive normal metal with weak diamagnetism. For 
T cz < T < T C 2 we have a vortex line liquid, with strong diamagnetism but still with resistive 
behavior in all directions. For T m < T < T cz we have a vortex line liquid with supercon- 
ducting coherence parallel to B. For T < T m we have an Abrikosov vortex line lattice. For 
T < T c ± the vortex line lattice is pinned. 

If we fit the lowest five data points (those for 77 < rj cr ) to a power law, we find T c ± ~ 
rr°- 88±0 - 09 , T cz ~ r7-°- 98±a05 , and T m ~ r7-o.66io.or_ The resultg for T ^ and ^ &re [n gQod 

agreement with our dimensional argument that characteristic temperatures at small rj should 
scale as T ~ t]' 1 . The agreement of T m with this form is much poorer. Whether this reflects 
the inclusion of too large values of 77 in the fit, or whether it reflects a poor determination of 
T m due to finite size effects or incomplete equilibration, remains unclear. At large rj, all three 
lines approach the constant value T C 2D , which we have found from independent simulations 
to be the melting temperature for an isolated two dimensional plane.0 

We see that T c2 for 77 > 77 cr becomes independent of 77, and is located at the same 
temperature as the specific heat peak in the ordinary {B = 0) 2D XY modellil (which lies 



about 10% above the 2D XY Kosterlitz-Thouless transition at Tk_t/J± — 0.9). Thus, at 
these high temperatures, rj CT does indeed mark the dimensional cross-over where our three 
dimensional system is behaving as effectively decoupled 2D layers; the cross-over T C 2 is due 
to the proliferation of vortex-antivortex pairs within these decoupled layers. However, at 
lower temperatures, we see no dramatic change in behavior for T c j_, T m , and T cz as Tj CT is 
crossed. Layers remain coupled, and in particular, while T c ±, T m , and T cz appear to merge 
as r] increases, we continue to find T c j_ < T m < T cz for all r\ > r\ cr studied. 

Note that in the limit of weak anisotropy, r\ — > 1, T cz and T C 2 become close, as was 
observed in earlier isotropic simulations.0 However, once the anisotropy r\ increases, T cz falls 
well below T c2 . The transition at T cz is thus clearly distinct from any mean-field- like cross- 
over phenomena. We will discuss this point in greater detail in the following section. Using 
the analogy between increasing 77 and increasing B as discussed in Section 0, the increase in 



the width of the vortex line liquid region in Fig. 10, as anisotropy increases, is in agreement 



with the general experimental features discussed in the Introduction. 

B. Vortex Line Fluctuations 

We now discuss several measures of the vortex line fluctuations in our model, in order 
to try and clarify the nature of the phenomena at T cz and T c2 . 

The first quantity we consider is A£^, defined as the the total number of vortex line 
segments due to fluctuations in direction ft, normalized by the total number of field induced 
"pancake" vortices in the xy planes, fL 2 ± L z . Note that in computing Ai^, line segments 
are added without regard to the sign of their direction; oppositely oriented segments do not 
cancel out. In Figs.|TT|a and |TT| fo we show our results for A£j_ = l(A£ x + A£ y ) and A£ z for 



the two cases of rj 2 = 10 and rj 2 = 50 respectively. We see that in both cases, A£ z is at 
least two orders of magnitude smaller than A£± in the vicinity of T cz and below. Thus, only 
transverse vortex fluctuations appear to be important at the phase transitions. Only at the 
higher cross-over T c2 , does A£ z start to become comparable to A£±. This is consistent with 
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our interpretation of T C 2 as the temperature at which vortex-anitvortex pairs start to enter 
the xy planes. 

One possible explanation for the transition at T cz has been proposed by Nelson,i in terms 
of the entanglement of vortex lines. If we assume, as in Nelson's picture, that the transverse 
fluctuation of a vortex line in the liquid phase is like that of a random walk, then since A£±_ 
is the net transverse fluctuation per pancake vortex, the total transverse deflection of a line 
in traveling down the length of the system will be u = y/L^A£±. Geometric entanglement^ 
of lines should occur when u ~ a v , or when A£± ~ a^j\fT7 z = fL z . For our system 
with / = 1/15, this criterion gives entanglement at values of A£± = 1.0, 0.71, and 0.35 
for thicknesses L z = 15, 30, and 120 respectively. Noting that A£± shows no apparent 
dependence on L z near the transitions, we would conclude that geometric entanglement 
takes place noticeably below T cz for systems of thickness L z > 15. It is interesting to note 
that in both cases Fig.|TT]a and |ll|fc, T cz appears to coincide with the point where A£± ~ 1. 
We have similarly observed this to be true at other values of rj. However we have no 
explanation for this coincidence. 

The above argument assumed that all of the vortex line fluctuations consisted of trans- 
verse motions of the magnetic field induced vortex lines. However, there is additionally the 
possibility of forming thermally excited closed vortex rings, which for large enough t] and 
temperatures low compared to T c2 should tend to lie between two adjacent xy planes. We 
now describe our algorithm to trace out the paths of vortex lines, which will allow us to 
measure both the distribution of such closed rings, as well as the entanglement of the field 
induced lines. We start by searching the plaquettes for a penetrating vortex line segment. 
We then trace its path into and out of subsequent unit cells of the grid. Such a line can 
belong either to a field induced vortex line, or to a closed vortex ring. Tracing the line, we 
measure the net displacement parallel to z that is traveled before the line closes back upon 
itself. If we have a closed ring, this net displacement is zero, and we measure the perimeter 
of the ring p. If we have a field induced line, then because of our periodic boundary con- 
ditions parallel to z, this net displacement must be mL z with integer m = 1,2,..., fL\. 
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If m = 1 the line closes back upon itself upon traversing the length of the system L z . For 
m > 1, the line belongs to a group of m lines that are braided with each other. This is 



schematically illustrated in Fig. [12]. The distribution of values of m is a measure of how 
geometrically entangled the field induced lines are. With this procedure, we search through 
all plaquettes until all vortex line segments are found, and classified as belonging to either 
a ring of perimeter p, or an entangled braid of order m. The only complication in the above 
algorithm occurs when two or more vortex lines segments intersect, i.e. go in and out of 
the same unit cell of the grid. In this case we randomly choose which segment is connected 
to which. In practice this was achieved as follows. Once a line was traced into a unit cell, 
we searched the remaining five faces in a random order to see which face the line is leaving 
through. Once we find a line leaving, we take it to be the continuation of the line we are 
tracing. 

In Figs.|13|a and |13fo we show our results for the distribution q(p) of the number of closed 
rings of perimeter p per unit volume L Z L 2 ± , for the two cases rj 2 = 10 and rj 2 = 50 (we show 
results for cooling; no significant hysteresis was observed comparing heating and cooling). 
Plotting the logarithm of q(p) versus 1/T, we see approximately straight lines at low T, 
indicating thermal activation. These lines have a change in slope in the vicinity of the 
melting T m , which is mild for rj 2 = 10, but more pronounced for rj 2 = 50. At the higher 
T c2 the curves saturate. We believe this is consistent with the interpretation of T c2 as the 
cross-over temperature at which, upon heating, vorticity explodes throughout the system, 
and superconducting order is lost on even small length scales. The saturation of q(p) occurs 
because, above T c2 , the distribution of ring sizes is governed more by the statistics of random 
intersections among the lines, rather than by energetics. Note that in Ref. [7|, where only the 
isotropic case was studied, we incorrectly associated this explosion of vorticity, as indicated 
by the saturation of q(p), with the transition at T cz . From the phase diagram of Fig. [10], we 
now see that this mistake was due to the proximity of T C 2 and T cz which occurs only in the 
isotropic limit. For anisotropic systems, T cz drops below T c2 and lies in the region where 
q(p) is still governed by thermal activation. 
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From the data of Figs.O and O, we can now compare how much of the vortex line 



fluctuations is contained in the wandering of the field induced vortex lines, versus how much 



is contained in the thermally excited vortex rings. We show in Figs.|T^a and |TJ6 the total 
length of all vortex line fluctuations (per pancake vortex), A£ to t = 2A£± + A£ z , and the 
total length contained in closed vortex rings (per pancake vortex), A4i ng = f^YlpPlip), 
for the cases rj 2 = 10 and r] 2 = 50 respectively. We see that at low temperatures, rings 
constitute a negligible fraction of the total vortex line fluctuations. At the transition T cz 
they are only about 3.5% of the fluctuations for rj 2 = 10, and 1.5% for rj 2 = 50. It thus seems 
that as anisotropy increases, the importance of disconnected thermally excited vortex rings 
decreases. It is worth noting however that ring excitations which connect to field induced 
lines, may still play a role in determining the wandering of the magnetic field induced lines 



between planes. An example is illustrated in Fig. 15. Such "connected" ring excitations, 



which lie between xy planes, are degrees of freedom distinct from the pancake vortices 



which lie within the planes. In Figs. |16p and |16p , we show snapshot views of vortex line 
configurations, at various temperatures, for the cases rj 2 = 10 and rj 2 = 50, for a system 
of size 15 3 . The bottom row in each figure is a view of all vortex line segments that lie 
between a typical pair of adjacent xy planes. A "U" shaped segment in these bottom row 
views indicates a connected ring excitation. We see that they are present in the system at 
virtually all temperatures shown, however it remains unclear how, if at all, they correlate 
with the transitions. 

We turn now to consider the fluctuations of the magnetic field induced vortex lines. 
As mentioned in connection with our line tracing algorithm, we can classify each such line 
as belong to an entangled braid of m lines, m = 1, . . . , fL\ (see Fig. 0). We denote by 
n(m) the number of lines which participate in a braid of order m, and R = n(l)/ fL\ is 
the fraction of unentangled lines. In Figs.[T7|a and [TT| fe we plot R versus T for the cases 
rj 2 = 10 and rj 2 = 50 respectively. We see that upon heating, R — 1 up to T m , throughout 
the vortex line lattice phase. Above T m , R starts to drop, tending to saturate to its high 
T limit around T cz . Upon cooling, R starts to increase at T cz , and in most cases reaches 

13 



a completely disentangled configuration with R = 1 at T m . These features were found at 
all values of rj studied. For the thickest sample of L z = 120 at rf = 10, however, we cool 
into a glassy entangled state with R ~ 0.47 frozen below T c ±, as has been seen in previous 
isotropic simulations.^ From Figs. [17] it seems clear that the transition at T cz is related to 



the braiding of lines. 

To see this another way, we plot in Figs. [l8|a and |l8l&, for rf = 10 and rf = 50 respectively, 



the braid distribution n(m) versus m, for several different temperatures near T cz . We use 
our data for systems of size 30 3 , which have 60 lines. We see that for T < T cz , n(m) is 
strongly peaked at small m, decaying rapidly as m increases. However, as T cz is approached, 
the peak at small m decreases and the distribution n(m) becomes flat and equal to unity for 
an increasingly wide range of intermediate m. The transition at T cz therefore seems to be 
associated with braids involving a macroscopically large number of lines. When n(m) = 1 
for all m, it indicates that a line i which starts out at ri±(z = 0) in the xy plane at z = 0, 
is equally likely to match onto the starting position of any other line j, after traveling down 
the thickness of the system, i.e. Vi±(z = L z ) = Tj±(z = 0) is equally likely for any pair i 
and j (see Fig. Yl). We may speculate that precisely this condition is achieved at T cz in the 



limit of large system sizes, and long simulation times. 

An intriguing question concerning behavior in the vortex line liquid is how easily lines can 
cut through each other. This has important consequences for line diffusion. If lines cannot 
cut, they can be effectively pinned by their mutual entanglement s.i@ For our system in 
particular, with periodic boundary conditions parallel to the magnetic field, the degree of 
entanglement can only change due to the cutting and reconnecting of lines. To investigate 
this we have computed the average number of vortex line intersections, N c , present in any 
instantaneous configuration of the system. An intersection is defined when two vortex lines 
enter and leave the same unit cell of the grid, and corresponds to vortex lines with overlapping 
cores. Once two lines intersect, they are free to cut through each other, or even to detach 
and reconnect different ingoing and outgoing segments. We define the "cutting length" 
£ c = fL\L z /N c as the average distance (in units of d) along z between cuts of the magnetic 
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field induced vortex lines. £ c gives a crude measure of the average length over which a vortex 
line remains a well defined string, or equivalently a measure of the number of planes which 
remain correlated. In Fig.[19| we plot £ c versus T for the two cases rf = 10 and rj 2 = 50. 
In the phase diagram of Fig. |1(] we show contours of constant £ c = 2, 4, 6, and 10. We see 
that planes are essentially uncorrelated at temperatures above the cross-over T c2 . However 
correlations grow and get large as one cools below T c2 towards T cz . The picture presented 



by the contours of £ c in Fig.|l0|, combined with the behavior of R in Figs.[17| and n(m) in 



Figs. 0, is that intersections start to freeze out below T C 2, with lines becoming well defined 
on longer and longer length scales. This presumably will affect the time scales on which 
lines are able to diffuse about, with a a corresponding signature to be expected in dynamic 
phenomena. However the equilibrium degree of entanglement, as measured by R and n(m) 
remains largely unchanged down to T cz . Below T cz , the behavior of R and n(m) shows 
that the lines start to disentangle, yet cutting is still frequent enough to change the degree 
of entanglement for all temperatures down to T m . Below T m the lines remain either in a 
disentangled lattice phase, or a metastable state with a frozen degree of entanglement. 



IV. DISCUSSION 

We have computed the phase diagram of a fluctuating type II superconductor in the 
anisotropy-temperature plane. Our results are consistent with general experimental ob- 
servations, that vortex lattice melting occurs well below the cross-over T C 2 associated with 
the formation of local superconducting order, and that the width of this region increases 
with increasing magnetic field (anisotropy). As was found in earlier isotropic simulations, 
we continue to find as anisotropy increases that there exists a finite region of the vortex 
line liquid which possesses superconducting behavior parallel to the applied magnetic field. 
Such a region has not been seen in other simulations^ which have used a higher vortex line 
density, / = 1/5, 1/6, and 1/8. The transition at T cz from the superconducting line liquid 
to the normal line liquid appears to be associated with the braiding of a macroscopically 
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large number of the field induced vortex lines, in qualitative agreement with Nelson's picture 
based on an analogy with the superfluid transition of two dimensional bosons.! However, 
unlike Nelson's picture, T cz does not appear to decrease with increasing L z . Such a pos- 
sibility has been proposed by Feigel'man and co-workers^ in terms of an analogy to 2D 
bosons with long range interactions. Recently, Tesanovic has proposedil a mechanism for 
such a transition in terms of a vortex loop unbinding analog of the transition at B = 0. Our 
numerical results indicate that disconnected thermally induced closed vortex rings do not 
appear to be playing any significant role at T C2 , once the anisotropy has increased enough 
that T cz is significantly below T C 2. However it remains unclear whether or not vortex rings 



between planes, which are connected to the field induced lines (see Figs. [15] and [16]), are 
important degrees of freedom. 

We have studied behavior as the anisotropy is increased beyond the "3D-2D" cross- 
over, which simple dimensional analysis gave as r]^ T f = 1. We found that this did indeed 
mark the cross-over from ^-dependent (3D) behavior to //-independent (2D) behavior, at 
the high cross-over temperature T c2 . i] CT also corresponds to the region where there is the 
widest relative width for the floating vortex line lattice, T c ± < T < T m . However we see no 
qualitative changes in critical behavior as r] cr is crossed. As r\ increases above t] CT , T c ±, T m , 
and T cz all approach each other, but we continue to find T c ± < T m < T cz . 

Glazman and Koshelev,! by considering the effect of elastic distortions of the vortex line 
lattice on interplanar phase fluctuations, have argued that T cz should decrease below T m , 
with a dependence T cz ~ rj^ 1 , when the anisotropy increases above a value ~ 10?7 cr . Frey et 
al.@ have argued that when rj 3> r] cr , the proliferation of vortex lattice defects will create a 
"supersolid" phase, leading to T cz ~ 1/lnry, which again falls below T m for large enough rj. 
Similar results were earlier proposed by Feigel'man et al.0 We were unable to equilibrate our 
system at such high anisotropies so as to more thoroughly check these predictions (very large 
values of r\ require relatively large values of L±, so that the total interplanar coupling energy 
remains large compared to T). Recently, simulations^ by Nguyen et al., of a vortex line 
model with a finite value of X± < a v , find evidence for T cz < T m for moderate anisotropies. 
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They find T cz ~ r]~ 2 , scaling with the coupling between planes, in contrast to the above 
two theoretical predictions. They attribute the effect as due to the proliferation of closed 
vortex rings lying between planes, in contrast with our own findings that such rings do 
not exist in significant numbers. The point where their T cz crosses below T m occurs at 
the value of anisotropy where magnetic coupling between planes starts to dominate over 
Josephson coupling, r\ > Aj_/£j_. Their results thus lie outside the range of validity of our 
infinite Aj_ model (see Eq.(||)). Simulations by Sasik and Stroud,0 on an anisotropic model 
in the "Lowest Landau Level" approximation (also a Ai -> oo approximation), always find 
T cz = T m for all anisotropics studied. Thus the possibility that parallel coherence can vanish 
at a lower temperature than melting, and if so the nature of the mechanism responsible, 
remain issues yet to be clarified numerically. 

Finally, we have shown that the effective length over which vortex lines can be considered 
well defined connected objects, as measured by the distance between intersections £ c , steadily 
increases once one cools below T c2 into the vortex line liquid. This is consistent with the 
analysis of "non-local" conductivity in flux transformer experiments, which indicate that 
correlations parallel to B start to grow right from the onset of strong diamagnetism.il 
Nevertheless, we find that vortex line cutting remains sufficiently easy over most of the vortex 
line liquid region. This is illustrated by the absence of any hysteresis in the measurement 
of our entanglement parameter R, for most of the temperature range T m < T (recall that 
for our periodic boundary conditions, R can only change value due to the cutting and 
reconnection of lines). Only for rf = 10 and our thickest L z = 120, did we find freezing into 
a non-equilibrium state with finite entanglement, that sets in near T m , and saturates below 
T c ±. In earlier isotropic simulationsB, this freezing out of equilibrium was found to occur at 
a higher temperature, a little below T cz . 
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FIGURES 

FIG. 1. Ground state locations of vortex lines in the xy plane for line density / = 1/15 on a 
cubic grid. 

FIG. 2. Helicity moduli Yj_ and T z versus temperature T for anisotropy rj 2 = 10 and vortex 
line density / = 1/15. Heating and cooling data for three different system sizes are shown, along 
with representative errors. 

FIG. 3. Structure function <S(kj_) for if = 10 and / = 1/15 for system size 30 3 , upon heating. 
The cross-over from Bragg peaks to liquid like rings occurs at T m /Jj_ ~ 0.43 

FIG. 4. Bragg peak heights A5(Ki)/5o for rj 2 = 10, / = 1/15, and different system sizes. 

FIG. 5. Specific heat C versus T for rj 2 = 10, / = 1/15, and various system sizes. The high 
temperature peak locates the cross-over T C 2. 

FIG. 6. Helicity moduli Yj_ and T 2 versus temperature T for anisotropy i] 2 = 50 and vortex 
line density / = 1/15. Heating and cooling for two different system sizes are shown, along with 
representative error bars. 

FIG. 7. Structure function 5(kj_) for rj 2 = 50 and / = 1/15 for system size 30 3 , upon heating. 
The cross-over from Bragg peaks to liquid like rings occurs at T m /J± ~ 0.21 

FIG. 8. Bragg peak heights AS(Ki)/Sq for rj 2 = 50, / = 1/15, and different system sizes. 

FIG. 9. Specific heat C versus T for rj 2 = 50, / = 1/15, and various system sizes. The high 
temperature peak locates the cross-over A lower temperature peak corresponds to T c ± 

FIG. 10. Phase diagram in the anisotropy-temperature plane for vortex line density / = 1/15. 
£ c is measured in units of d. T C 2 locates the peak in specific heat. 
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FIG. 11. Excess vortex line length due to fluctuations in transverse, A£j_, and parallel, A£ z , 
directions, a) is for rf = 10, b) is for rf = 50. Dashed lines are guides to the eye only. 

FIG. 12. Schematic example of the possible reconnections of field induced vortex lines, under 
application of the periodic boundary condition in the z direction. Solid, dashed, and dotted lines 
are used to distinguish the different lines within a particular braid. 

FIG. 13. Number of closed vortex rings q(p) of perimeter p, per unit volume, a) is for rf = 10, 
b) is for rf =50. 

FIG. 14. Total length of all vortex line fluctuations, A^ot, and total length of lines in closed 
vortex rings, A£ r ; ng , per number of pancake vortices, a) is for rf = 10, b) is for rf = 50. 

FIG. 15. Schematic example of how "connected" vortex rings between planes contribute to the 
wandering of field induced lines between planes. 

FIG. 16. Snapshot views of vortex line configurations at various temperatures for a) rf = 10, 
and b) rf = 50, for 15 3 size system. For each case we show perspective views from the side (top 
row), and looking straight down along the applied field (middle row). Solid dots indicate a point of 
intersection between two vortex line segments. The bottom row is a view of all vortex line segments 
that lie between a typical pair of adjacent xy planes. 

FIG. 17. Fraction of unentangled lines R, vs. T, for various system sizes, a) is for rf = 10, b) 
is for rf =50. 

FIG. 18. Braid distribution n{m) vs. m for several temperatures near T cz , for system size 30 3 . 
a) is for rf = 10, b) is for rf = 50. 

FIG. 19. Cutting length £ c vs. T for rf = 10, and 50, for various system sizes. 
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